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ABSTRACT 

We propose a unified description of tfie galactic Dark Matter and various 
uniform scalar fields for the inflation and cosmological constant by a single field. 
The two types of effects could originate from a fluid of both spatially and tem- 
porally varying Vacuum Energy if the vacuum has an uneven pressure caused by 
a photon-like vector field (of perhaps an unstable massive boson). We propose 
a most general Lagrangian with a Non-uniform Cosmological Constant for this 
vacuum fluid (dubbed as a Nu-Lambda fluid), working within the framework 
of Einsteinian gravity. This theory includes a continuous spectrum of plausible 
dark energy theories and gravity theories, e.g., inflation, quintessence, k-essence, 
f(R), Generalized Einstein- Aether f(K), MOND, TeVeS, BSTV etc. theories. It 
also suggests new models such as a certain f(K+R) model, which suggests in- 
triguing corrections to MOND depending of redshift and density. Some specific 
constructions of the Nu-Lambda fluid (e.g., Zhao's V-A model) closely resemble 
the ACDM cosmology on large scale, but fit galaxy rotation curves as good as 
MOND. Perturbed Einstein Equations in a simple /(-ft'4) model are solvable and 
show effects of a Dark Matter fluid coupled and Dark Energy. Incorporating the 
perturbation equations here into standard simulations for cosmological structure 
growth offers a chance to falsify examples of the Nu-Lambda theories. 

Subject headings: Dark Matter; Cosmology; Gravitation 



1. Introducing a framework for vector fields 

General Relativity (GR) is actually a special case and minimal construction of a range 
of theories describing the metric of a plausible universe. While completely adequate on small 
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scales, GR by itself predicts a missing mass and missing energy compared to astronomical 
observations of the metric of the universe on the scale of Kpc to Gpc (e.g., Spergel et al. 
2007). While the missing mass is arguably explained by Dark Matter (DM) particle fields 
in supersymmetry particle physics, the missing energy almost certainly cannot be explained 
unless the present universe is immersed in an exotic Dark Energy (DE) field (White 2007). 
Since both the effects of DM and DE occur when the gravity is weak, one wonders if the 
underlying fields are tracking the metric field of the gravity (Zhao 2006, 2007, 2008). 

Quantum gravity and string theory often predict a non-trivial coupling of some vector 
field, which violates CPT symmetry satisfied by standard physics (Kostclccky & Samuel 
1989, Kostelecky 2004). It has been considered by Will & Nordvedt (1972) that a vector field 
can be coupled to the space-time metric. This creates a "preferred frame" in gravitational 
physics. A global violation is undesirable, but a local violation is allowed. A four time-like 
vector-field with a non-vanishing time component would select a preferred direction at a 
given space-time coordinate. It is an aether-like fluid present everywhere, somewhat like a 
dark energy with some preferred direction. If such a vector coupling to matter is zero or 
small, then it can evade current experimental detection (e.g. the CPT violation experiments 
in Princeton). There have been an increase of interests about such vectors in recent years, 
especially works by Kostelecky and coworkers (1989, 2004), Foster & Jacobson (2006) and co- 
workers, Lim (2005) and co-workers, Bekenstein (2004), Sanders (2005), and Zlosnik, Ferreira 
et al. (2006, 2007) and co-workers. E.g., Foster & Jacobson (2006) noted that a solar system 
immersed in a unit time- like vector field (called Einstein- Aether, or AE) of small enough 
mass coupling to the metric is apparently consistent with current measurements of PPN 
parameters. Carroll & Lim (2004) noted that such a field can have effects on the Hubble 
expansion. Inspired by these ideas, several workers especially Bekenstein (2004), Sanders 
(2004) and Zlosnik et al. (2007) proposed to extend the application to galaxy scale to use it 
to explain missing matter (Dark Matter). Many have constrained the theory using empirical 
astronomical data (Famaey & Binney 2005, Zhao & Famaey 2006, Zhao 2006), including 
gravitational lensing (Zhao et al. 2006, Chen & Zhao 2006, Chen 2008, Angus et al. 2007). 
Most recently Zhao (2007) found a simple Lagrangian within these frameworks to give rise 
to the DM-DE effects of the right amplitude, offering a possible explanation of coincidence 
of DE scale and DM scale in ACDM cosmology. The model is dubbed Vector-for-Lambda, 
or Vacuum-for-Lambda, or VA. 

Here we propose a very general Langragian of a non-uniform vector field in the vacuum. 
It is effectively a Dark Energy fluid in the GR language. This fluid with a Non- uniform 
Lambda is dubbed a Nu-Lambda fluid. We show its relation with existing theories. We 
isolate a simple case, and give the full field equations. Most importantly we derive the 
equations governing perturbation growth in the FRW universe. 
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To build a covariant theory, one starts with the Einstein- Hilbert action used for GR 



S 





R 







+ Sm [l] 



where the hght speed is c = 1, and Gjq is the gravitational constant, and g is the determinant 
of the metric gap] The signature taken here is (—,+,+,+); we do not distinguish Roman 
ahcd and greek ajd-^d for four indices. R is the Ricci scalar, describing the curvature of 
space-time. Sm is the matter action that describes the matter distribution. Variation of this 
action with respect to the metric gives the Einstein equations: 

^ = SvrGT™/"^'- (2) 



where Ga/3 = Rap — ^QapR is the Einstein tensor, and T^"**"^^: the stress-energy tensor of 
matter defined by: 6Sm = ~| / d'^Xyf—gT^i, (x) 6g^'^ (x). This tensor describes the matter 
distribution. In the following we will add new fields and terms to the Einstein-Hilbert action. 



1.1. Vacuum vector field with a dynamic norm 

Denote a vector (Bosonic) field in vacuum by Z°', which has generally a variable or 
dynamic norm. The Lagrangian density of many vector theories can then be casted in the 
general form0 



L = il + fo{ip))R + fi^,K,J) (3) 

= Z-Z, (4) 

K = KfaVaZ'V.Z" (5) 

J = J^^VaZ'VcZ'', (6) 



where cp'^ is essentially a scalar field made from the norm of the vector field Z'^ without 
introducing new degree of dynamical freedom, the coefficients K^j^ and J^J' can be lengthy 
functions of and the metric gab with appropriate combinations of upper and lower index. 
For example, the Lagrangian of the VA model (Zhao 2007) is of the form 

L = R + fK{K) + fj{J) (7) 



^It is optional to add a new term f4{ip)Z'^Z^Rab, which is related to f{(p, K, J) via a full derivative of no 
effects. 
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where K = {Z''VaZ''){Z''VcZ'^), J = 5^52{VaZ''){VcZ'^). To recover scalar-tensor theories 
or f{R) theories for Dark Energy, we set coefficients of K and J to zero, so end up with a 
Lagrangian 

L = R + fo{^)R + f{^), (8) 

where the vector field Z°- has collapsed into its norm, the scalar (p. All these thoeries involve 
only four degrees of dynamical freedom at maximum. 

As a specific illustration where the coupling coefficients are the simplest, the Lagrangian 
scalar density can be 

8 oo 

L{Z,g) = R + Lou + L3 + f^U + a^U + J] aiU. (9) 

i=5 4=9 

For simplicity consider setting all coefficients = for z = 9, ...00, we have 

-^012 = flO + O-lV^ + Ci2V'^-, = Z'^Z, 

L3 = asip^R 
Z/4 = Z'^ Z^ R^fs- 

and 

L7 

where the coefficients Oj = const. 

This Lagrangian density can be simplified further; note the L4 term is related to Lg and 
L7 by a total divergence. 

M4 = [Va(/4PF") - ly^Va^] - /4L6 - /4L7 (17) 

iy« = [Z'^VbZ^ - Z^VbZ") , (18) 



(10) 

(11) 
(12) 



(V,Z^)(V"ZO 
{V^Z^f 
(VaZ^)(V^Z") 
{ZNpZ^){Z^V^Z^ 



(13) 
(14) 
(15) 
(16) 



^The Lagrangian of the BSTV theory of Sanders (2005) can be casted into a simihar expression but with 
5 degrees of freedom in the physical frame, with L = R + d{(p)g'^''\/aq'^bq + h{q, (p)K — f{q, ip)J + 2V{q, ip), 
where g is a new dynamical scalar field. In the slow-roll approximation, we can neglect the dynamical term by 
setting d{ip) ~ 0, eliminate q all together by minimizing the action wrt q, hence rewrite L = R + f{ip, K, J). 
Essentially the function / is a slow-varying scalar in BSTV with its five degrees of dynamical freedom. 
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here we can drop the term proportional to V/4, which is zero if /4 is a(/7-independent constant. 
The total divergence term, when integrated over volume, can be dropped in the total action 
S, because the term becomes a surface integration over the boundary according to the Stokes 
theorem. We can therefore choose not to consider the term L4, absorbing its contribution in 
the Lq and L7 terms. 

1.2. A general Lagrangian of a scalar field plus a unit- vector 

So far Z" is a vector, not required to be unit-norm. It is, however, easier to work with 
unit vector. Now we decompose the vector Z°- into a scalar field (p representing its norm 

= ^JE", (19) 

pis a unit time- like vector Ma — Z^j Basically 

= ^\ M'^Ma = 1. (20) 

Note the sign convention for the unit vector. 

The time-like unit vector guarantees Lorentz invariance to be broken locally, so that 
it will always have a non- vanishing time-like component. The additional constraint can be 
enforced using a non-dynamic Lagrange multiplier L* . 

The co-variant derivative Z-terms have the following correspondences to co-variant 
derivatives of the M field and the scalar (/? field: 

{Z''VaZ,){Z'VbZ')^Ls = ^'"M^M'^Va'pVb'P + 'P^K^ (21) 
where the KiS are defined as 

= M''M''g^VaM''VbM'' = M'^VaM.M'VbM''. (22) 
Inspired by the above specific case and redefining the coefficients, we propose a very 
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general Lagrangian for a dynamical scalar field (p coupled with a unit-norm vector field iE", 

4 

£((^,iE) = [l + co{^)]R + 2V{y^) + J2c^{^)K^ + {JE-JE,-l)L* 

1=1 

+ [dsiip)M''VaM'' + d^{if)M'VaM''] Vb^ (23) 

where L* is the Lagrange multipler, R is the Ricci scalar, and c's and d's are now treated as 
general functions of the scalar field ip and these terms are some kind of dynamical dark energy; 
(y9 is a singlet (a real number). It can be turned non-dynamical if di = d2 = d^ = d4^ = 0. 
Models with d2 = d^ = d^ = are simple quintessence models. 



1.3. Nu-Lambda fluid vs. other theories 

Our general Lagrangian L{ip, JE) for the Nu-Lambda fiuid includes the Lagrangian 
L{Z,g) in eq. ([9]) as a special case. In fact all coefficients in L{Z,g) could be gener- 
alized to functions of (/?. To see this note that f4{ip)Z°'Z^Rab in eq. f[T7|) would contain a 
term proportional to 

W'Vah = flW'^Va'P (24) 

and W"' can be casted into and IE, and the end result are terms all included in L{ip,JE). 
Likewise, any new terms ajLj can be collapsed into and IE representation. For example, 
it was proposed (Ferreira et al. 2007) to set 02 = as = to eliminate the a^Rp"^ coupling 
term, but include four new terms OgLg -f- aio-Z^io + + 012-^12 in eq. ([2]), would not lead 

to new terms in our L{ip, M). To see that how these new terms are absorbed, we note that 

Z'Z.WaZ'WbZ'' 
Z^Z.VaZ^VbZ' 
Z,ZdWaZ''W''Z'' 
Z^^Z'Z.ZdVaZ'WbZ'' 

Note the right hand sides are all already included in our Lagrangian eq. (l23l) . with a specific 
assignment of our functions Ci, 2,3,4, c^i, 2,3,4- § 



= Lg = p)'^IE''M'S/a'pS/bV + ^'^M'S/bIE°-Va'P 

= Lio = ^''M^'IE'VapVb^ + p^M'^Va^V^'Mb 

= Lu = pi^g'^^VaP^bP 

= L^2 = P^IE''IE'WapyhP. 



•^Oiir Langrangian for is a special case if we choose cq — 03(^9^, — o-s.ej'p'^i C4 = as^p^, rfi = 

ay + 011.^2,^2 = 05 + fle + (as + 09 + aio)^?^ + a^ip'^, da, 4 = Sas^g-^ + ag^iQip^. 
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Bekestein's (2004) TeVeS Lagranian can also be casted (Zlosnik et al. 2006) into that 
of a pure non-unit norm vector field in the physical metric with a Lagrangian 

L = R + fj{J) + fK{K), (25) 

where fj{J) = J and the functional form of fxiK) is determined the MOND interpolation 
function. The variables K and J are made of terms 

12 12 

K = Y,U^)Li, J = ^Ji(<^)^i, (26) 

where hi and ji are functions of the norm ip, and Lj are the 8 different combinations of the 
kinetic terms of the non-unit norm vector Z", which then reduces to our Lagrangian L[ip, M) 
using scalar field and the unit vector. 

In fact we claim that our Lagrangian I/((^, M) is general enough to include several models 
in the literature as its special cases: 

1. GR: GR corresponds to our model with 6^ = 0, Oj = except that 02 and ao- 

2. Scalar-tensor gravity: when Ci 2,3,4, ^2,3,4 = and rfi = 1 it reduces to the scalar-tensor 
gravity; if furthermore Cq = it becomes the standard scalar field theory for infiation 
and quintessence etc. 

3. f{Ft) gravity model: when Ci^2,3,4, ^^1,2,3,4 = we could vary the action with respect to </? 
{nondynamical now) to have = ^ CQ{(p)R + V'{ip) = where a prime means d/d^p. 
Solving this equation we get ip{R) so that the action becomes that for the R -\- f{R) 
theory. 

4. Einstein- Aether model and f{JC) model: set Cq, (^1,2,3,4 = 0. When Ci^2,3,4 = ^"1^2,3,4 are 
constants and V{(p) — one obtains Jacobson's linear ^-theory (Jacobson & Mattingly 
2001) with a. More generally assuming £1,2,3,4 = Ci,2,3,4</', we can again vary the action 
wrt ip and obtain Yli CiKi + V'{ip) = 0. We can solve p as function of /C = CiKi, and 
write the Lagrangian as /(/C) (Zlosnik, Ferreira, Starkman 2007, Zhao 2007). 

5. V-A model: This requires a non-dynamical scalar doublet (A^, Aj) (Zhao 2007). We 
set (ii, 2,3,4 = co,i,3 = 0, and C2, C4 are functions of the two independent components 
(Ak, Aj) of the doublet respectively. 

6. TeVeS: This requires the scalar field to be a doublet {(p^pL), that is, it requires two 
scalar fields. The field is non-dynamical (cf. Bekenstein 2004, Zlosnik, Ferreira, 
Starkman 2006). The expressions for our functions are lengthy. 
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7. BSTV: This again requires a scalar doublet (</?, q) (Sander 2005). But we set di + (^2 = 
h{q) and = /(?), and ^3 = ^4 = Co,i,2,3,4 = 0. 

There are various other special cases, not in the literature, e.g., If we set (ii, 2,3,4 = 0, 
ci,2,3,4 = <^C'i,2,3,4, where C's are constants, and co(0) 7^ 0, we can vary wrt tp to get an 
equation of motion, (/? = (p{R,K), and eliminate the non-dynamical (f and its potential 
y (</?), and then cast the Lagrangian as L — R + F{R, K) models. This kind of models are 
simpler than TeVeS, since the unit vector has only 3 degrees of dynamical freedom and there 
is no scalar freedom. 



2. The F{Q) models 

The dynamics of our general Lagrangian is very rich. To be specific, let's consider 
the simpler F{Q) where Q — coM'^R + JC models, where we redefine Cj as constants for 
i — 1,2,3,4, and the Cq term includes a dimensionless (linear) dependance on the Ricci 
scalar, and redefine variables so that the final total action is 

S = Sm + J d'x^g [R + M^F{Q) + L* (A"A« + 1)] (27) 

Q = cqM-^R + K, (28) 

where L* is the Lagrangian multiplier, and 

IC = M-'^K"^ V^A^V(iA% (29) 
i^"^ 7- = ciy"^^-ya + C25^5^ + C35^5? + c^A'^Afg^,. (30) 

Note that we have replaced the M field with A field using the opposite sign convention. 

A'^A^ = -I. (31) 

We will stick to A field (instead of M field) for the rest of the paper. 

In the case that Cq = 0, our action is similar to what was considered by Jacobson and 
co-workers, except for the non-linear F-function. Notice that dropping the terms in C2 and 
C4 and considering C3 = — ci, we find K°'l^ aA'^V pA'^ — ^FaaF""^, where Faa is the 
antisymmetric Maxwell tensor defined by F^a — Vq^o- — VaAa- This simplification was 
used by Jacobson and by Bekenstein in TeVeS. 

Models with 01,2,3,4 = and cq 7^ are F{R) theories. Models with cq = C4 = have 
been studied by Zlosnik et al. (2007) without giving the full equations. Here we expand on 
previous results. 
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2.1. Fields equations for F{Q) models 

Now we proceed to obtain the field equations for models, where Q = cqM~'^R + /C. 
What must be borne in mind when carrying out the variations is that the two dynamical 
degrees of freedom considered are the inverse metric g'^" and the contravariant vector field 
A^^. The contravariant vector is chosen (and not the covariant one) just because once one 
has chosen to variate the action w.r.t. g^'^, the result of this variation will be simpler seeing 
the form of K"'^ ja- because we have : 

= (32) 



bg 



where we used the fact that 



SwS'" = ^l^ hi« = -9w9-o'5»'" (33) 



The vector equation is obtained by varying the action w.r.t. A^: 



= ^ V„(F' J" - F'yp = 2L*Ap (35) 



5A 

which we define F' — ^ and 

• J"^ is a tensor current: J"^ = {K'^^^ +K^'^ ^a)^ = ^K'^^^V pA^ , due to the 
symmetry here in K. 

. = V,AW^>1«^^^^ = c^M-'^A^Vo^A^VpA'^. 

To get the Lagrange multiplier L*, we multiply the equation by and contract. Once 
L* is known, the equation (which has four components /9 = 0, 1, 2, 3) yields three constraint 
equations for the vector. Varying the action w.r.t. L* will give the constraint on the norm: 
A^A^ = -1 

For the variation of the action S = J d^x^—gL w.r.t. the contravariant metric, one 
must notice that = / d'^x^/^ (^P^" ~" \9(^I3^ where one uses the fact that: 6g — 

gg^^Sg^i, = —gg^vSg^^ , g being the determinant of the contravariant metric. 

The symmetry of K°'^ simphfies the equations: 

Sga/3 = WaP + F[Ya)3 + J ,—^^] (36) 
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with ^ 

W^f,^{ER^p + g^pWE-V^VpE), E= ^^^^\ (37) 

and 

i'., = V,A"V-,A^'-^^^. (38) 

The variation of the covariant derivative of the contravariant vector field requires to 
vary the Christoffel symbol (only) : 

5{V.A^) _ 5{d^A^ + Tl,A<^) _ 5{Ti;) 

SgaP 5gaP Sg^P V-^^J 

And we have 5 (F^^) = ^ {^a^gpr + ^ p^Qar — ^rSgcrp), so one eventually find 

F'J" r, ^-^J^ = -^V.(^'( J(„ Ap) - J\^Ap) - J^^p^A'^)) (40) 

dropping divergence terms which would once more contribute only by boundary terms. The 
brackets denote symmetrization, ie for instance J(ai3) = \ {Jap + J^a)- 



Putting these together and use 

SAi'A 



-A^Ap (41) 



5gOcl3 

we find 

= SttGT-/**- + f^p - W^p (42) 
and is the stress-energy tensor of the vector field 

Tad = 2^<^(-^V(a''^/3) - J^iaAfj) - J(^af})A'')) 

-:F'Y^ap) + \9apM-'T + L*A^Ap. (43) 



The above equations are actually true for all F(i?, K) models, if F' is interpreted as 
dF/dJC. It is interesting that the effect of the cqR term behaves partly as a rescahng of the 
gravitational constant by a factor (1 + S). Like in F(R) gravity the value of S oc Cq must be 
very small, to prevent the term Qab^'^E in the correction source term Wab to violate stringent 
constraints on small scales, e.g., the solar system. Unless stated otherwise, we will set this 
source to zero for simplicity. 



3. Full equations for perturbations of -F(/C) models in an expanding universe 

3.1. Metric, Matter, and Einstein Tensor 

With above equations of motion we consider a Friedman-Robertson- Walker (FRW here- 
after) perturbed metric such that 

ds^ = -(1 + 2e(t))dt^ + a{tf{l + 2eip){dx^ + dy'^ + dz^) (44) 

A static universe is a special case with all quantities independent of time and a{t) — 1. This 
fairly general metric is a weakly perturbed form of a homogeneous and spatially isotropic 
universe ds^ — —dt^ + a{ty{dx'^ + dy^ -\- dz^) by setting e = 0; the unperturbed metric is also 
spatially fiat here, i.e. with no curvature parameter. The potentials and ip are Newtonian 
gravitational potentials, which are generally non-identical. This form of perturbed metric 
neglects both tensor mode (gravitational wave) and vector mode perturbations. 

In the following, the equations are developed in orders of e, but e is not kept for a 
lighter expression. 

Matter: For matter fields, we can take 

which is the stress-tensor of a perfect fluid without any anisotropic stress, with a density 
p, a pressure P and with the fluid four- velocity satisfying g^j^u^u''' — —1. If we consider 
a non-relativistic fluid, hence neglecting the spatial components of u^, then in our metric 
U/j, — (—1 — 60, 0, 0, 0). We have also 

rpmatter ^ ^ 20) p 

rpmatter rpmatter n 

-'Oi ^ ■'■ij ~ ^ 

rpmatter ^ ^2 ^1 + 2^) P 



(46) 

(47) 
(48) 



Einstein tensor: Up to linear order in e we find 

. Goo = 3i/2 + QHdtij - ^ 

• Goi = 2{Hdicl) - dtdi^P) 

• G^^ = (a 2 +2a a) [-1 + 2(0 - i;)] + {dy ^ +d, 2) (0 + V') - 2a?dt ^ V + 2a a 9^(0 - 3^') 

• Gij = -didj{(j) + ip) for i ^ j 
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3.2. Vector field 

We take a homogenous and spatially isotropic universe for the background, so the vector 
field must, in the background, respect this isotropy for the modified Einstein equations to 
have solutions, so only the time component can be non zero. The constraint on the norm 
is QaisA'^A^ = —1 so in the background, we take A" = Sq and one can then expand it and 
write: 

A^ = g^^A"" = (-1, 0, 0, 0) + (-60, eA,, eAy, eA,). (49) 

The constraint on the total vector qqqA^A^ ~ — 1 and qqq — —{1 -\- 2e0) with the perturbed 
form of the metric fixes A^ 1 — ecf) and ~ fi'oo^'^ ~ — 1 — e0. 

We can also derive VA up to linear order in e, whose non-vanishing components are 
ViAo^-HAi, VoA - V^Ao = 9i0 + A, ViAj^a^ V' + ^(1 + 2V' -</>)] (^y + 9i>lj-(50) 

Following we carried out the calculations of the Einstein equations up to linear order 
analytically with this metric. Wherever the expressions become very lengthy, it is helpful 
to break the expressions into a non-spatial part without the Ai terms and a spatial part with 
Ai terms. Note also that any term F and its derivatives F' and F" contain implicitly an 
unperturbed part and a perturbed part. Finally we define the shorthands a = Ci + 3c2 + C3, 
and define df = didi, where i — 1, 2, 3 and di is the co-moving spatial derivatives and 
X — dfX is the synchronous cosmic time derivative. 



3.3. Kinetic scalair and F{]C) 



We decompose the kinetic scalar /C into the 0th, 1st and 2nd order terms. 



51C = K- 



M2 



M2 



-20 + 2//- V + 



3am 



diA, 



C4 - Ci 

a2M2 



where ... includes the 2nd order terms {-AH^jOtij + ^H^^ - 4:H(f)dtij + {dti^f)+^<t>dt<t> 
and other lengthy terms shown elsewhere (Halle 2007). All 2nd order terms are negligible 
for the linear perturbation calculations. However, the very first term (9i09j0) should 
be considered for static galaxies, where are in the non-linear regime. 

The terms F, F', diF'{K) are often involved in the Einstein Equation. We note they 
have different orders of magnitude, given by following 



5F = F-F 



M2 



F'SIC, 



(52) 
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5F' = F'-F 



f"5k:, 



diF'-O = F"{IC)dilC, 



6a 



where we have moved the 0th order terms to the LHS. In the special case where a = 0, we 
find SF', 6F, diF', 6{IC) are all zero up to the 2nd order. 



3.4. The Lagrange multiplier 

The vector equation gives the Lagrange multiplier L* — L*^ + L*^, where 



- *A 



3F' 
1 + 20 



(ci + C2 + C3) {H^ + 2Hi;) + C2( — + H(P- 6^^) 

a 



-F'{H + 3i:)-^{F'd,(P){53) 
a 



= +2a^F'd,A - ?>C2^d,{F'Ai) {aF'd.A^) - (F'A^ (54) 

a'* a"* a-^ V / 

where wc could drop a 2nd order term —2>C2-^ AAF' — ^diF'Af, note that both diF' and 
A are of 1st order. 



3.5. Perturbed equation of motion of vector field 



For the j-component of the EoM of the vector field, we have 



n C4 - Ci 
= dt 



aF'{dj(p + Aj 



+ 



C1 + C2 + C3 , 



F'^i^^Aj 



(djip - Hdj(l))F' + HdjF' + Ajdt{HF' 



(55) 



where we have dropped 2nd order terms involving product of {diF') and another first order 
quantities (^j^'jAj). This equation resembles equation Bl and B2 in appendix of Lim 
(2006). We consider only the scalar mode here Aj — djV for j = 1, 2, 3. 



3.6. The spatial off-diagonal terms 

The EE and vector field stress term with i ^ j satisfies up to linear order 

Cl + C3 



2a 



-dt[aF'{diAj + djAi)]. 



(56) 
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We can, as in the static case, identify the Newtonian potentials, + -0 = in the absence of 
the anisotropic stress T-^, i.e, in the (magnetic) case Ci + C3 — 0, e.g., GR is such a special 
case. 



3.7. The Oi Cross terms 

The Oi component of the stress tensor 

foj = ^^^dt (aF'{d,cj> + Aj)) + adt {F'H) Aj + A, (57) 

where A = +^^^di {F' {djAi — diAj)) is a curl-like term, and could be dropped in case of 
the scalar mode where Ai = diV; even for vector mode, one can drop the 2nd order term 
^^^diF' {djAi — diAj) safely. So we have the Ox-component EE 

+ 2Hdj(f) - 2dj'^ - {^^^^^ dt{aF\dj(t) + Aj)) = -SirGpuj - (58) 
for a non-relativistic matter fluid. 



3.8. The 00th Einstein Equation 



Replacing L*, we find the 00 component of the vector field stress-energy tensor Too 
satisfies 

^ C4 — Ci 



-00 



-d, 



a" 



{di(j) + Ai)F' + aF'H{3H + 6^ + 2a-^diAi) - 



^±^M^F (59) 



where we could drop a 2nd order term ^^-r^diF' A^. 



Thus we can write the 00 Einstein equation with T^"**^^ = (1 20) p as 



3(1-«F') (h\1 - 20) + 2H^)-'^^^diA~-^d1i: 



where we moved the F term to the RHS, and divided the factor (1 -|- 20). 



■^^a,(F'(a,0+A)) = U-kGp - ^ 



(60) 



3.9. Spatial Diagonal Equations 



The spatial diagonal terms satisfy, e.g.. 



7;, = a2(l + 2V^) 



1 - 



—a- 



-dt{F'a'd) + 



M^F 



+ aa 



-F'lp + F' (h4) - 6Hip - a^V) + A^ + Al(61) 
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where 



a Hdi{F'Ai) + -^{aF'diAi) 



A2 = 



Cl + C3 

3a 



dt [aF'{3d,A, - diAi)] (62) 



where the summation over i = 1, 2, 3 is imphcit, and we could drop a 2nd order term 
-aHdiF'Ai. 

Since for matter T™''**^'" = (1 + 2iIj) P, the modified pressure equation by adding the 
three spatial diagonal equations becomes 



where we have moved the vector field term to the LHS and divided the factor (1 + 2'0)a^ on 
both sides of Einstein's Equation. 



We have thus obtained the perturbations of the vector field stress-energy tensor and 
the Einstein equation for an vector field with a Lagrangian involving a general function of 
the kinetic term /C. As first check, wc recover the linear F — K, model of Lim (2006), and 
extend it to include a C4 term; this is given in Appendix. These perturbation equations are 
also consistent with Li et al. (2007), which use a very different formulation. As a summary 
of equations and further illustrations, let's consider some more special cases in the context 
of Dark Matter and Cosmological Constant. 

Two important quantities of later use are 



As it will be evident later on, this choice of F{1C) recovers MOND in present-day galaxies. 




4. Special Cases 





(64) 



4.1. F(/C4) models with ci = C2 = C3 = 



The perturbation equations become much simpler if we concentrate on models with a 
pure C4 term. By letting ci — C2 — — 0, we neglect all contributions of other kinematic 
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terms (one can set 0^ — 2 with no loss of generality). Up to the linear order the Lagrange 
multiplier 



L* = 0. 

We find t/j — —(f) from the spatial off-diagonal EE 

Gij = -9,9,(0 + ^) = f.^ = 0, for t ^ J. 



(65) 



(66) 



Collecting terms in the equation of motion and replacing — —0, and considering only 
scalar mode A^ — diV we find a time-independent quantity 



Xi = [Ma{tf]\E,, Ei 



Ma{t) 



84 + A 



(67) 



The quantity behaves as a time-independent co-moving gravitational force of an effective 
"dark matter" perturbation. It is the product of a dimensionless "electric" part by a 
dielectric " suspectibility" part A, and a changing scale of gravity Ma(t)^; note E^ is similar 
to the Electric field in a four potential (0, Ai, A2, A3) of the radiation field. 

The "suspectibility" or the dimensionless /C is an implicit function of the modulus 
of the dimensionless " electric" field E^ through the expression of 



\E\ 



/C 
2' 



A 



die' 



(68) 



where we set C4 = 2. 



The spatial variation of the time-dependent vector field A^ tracks the time-independent 
terms di{a(t)) and Xi by the EoM constraint 



a 



i - di{a4i) 
A 



(69) 



Time-wise, as the universe expands, /C, all approach 0, and A approaches a finite value, 
hence Aj — > 0. 



The stress tensors of the non-uniform dark energy fiuid are given by 

2dtX, 







= -A 



(70) 
(71) 
(72) 
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where A is a time-dependent effective "cosmological constant" term, and pDM,com is a time- 
independent effective co-moving "Dark Matter" density. 

The Einstein Equations for Gq?, Gq and become 

-dtdi(a(l)) = -SirGpUir^O, (73) 

a 

3H' + ^d^<l> = +A + SttG (^^^ + p) (74) 



a 



(1-20) = -A + SttCP. (75) 



where we assumed the source being the non-uniform dark energy stress tensor T plus a non- 
relativistic baryonic perfect fluid of pressure P and density p — {1 + 5)^^, where 5 is a 
growing overdensity, and Pcom is a constant co-moving mean density of baryons. 

In general the vector field = (—1 — 0, Ai, A2, A^) in F{lCi) models simply tracks the 
space-time metric perturbation and scale factor a(t), which tracks the dominant source, 
be it radiation or baryonic matter. Metric perturbation can be printed in the ^1,^2, A3 
fields even in the absence of baryonic matter. Note that the effects of the vector field is 
more complex than a change of gravitational constant of a baryon-radiation fluid, where 
Silk damping can erase perturbations; the vector field is not coupled to photons or baryons 
directly, hence its perturbations can be passed onto baryons after last scattering. The Tqo 
stress contains a DM-like source term, which decays with the redshift as fast as the baryonic 
density p = Pcomd'^, but keeping the effective DM-to-baryon contrast time-independent, i.e., 

- — P^^'"^'^ — independent of time, ^T^GpDM,com = diX^. (76) 

^ Pcom 

We can further introduce another parameter for the equation of state parameter defined 

by 



TO 



1 + 



A 

Clearly in case of early universe and CMB, is big, so 



(77) 



w = 0, (78) 

i.e., the equation of the state of the vector field is almost exactly Dark Matter like. This is 
important to understand why the vector field can replace Dark Matter (DM) in galaxies. We 
will show next that the F{K,^ model is essentially a non-uniform dark energy. A difference 
with real DM is that DM density perturbations can grow in co-moving coordinates while the 
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"dark" source term in F(/C4) corresponds to a static non-uniform density in co-moving 

coordinates. These effective DM and DE are actually coupled through the function JF, i.e., 
A is a function of pDM,com/ di^Y ■ E-g-j in the limit of weak field, K, = 2\E\'^ ^ 0, A — 1, 
A ~ Ao - M2/C/2 = Ao - {Xf + X| + Xi)a{t)-\ and pDM,coma-^ = ''"''t^Gai^'"''' ^ ^^ere 
{ki, k2, k^) is a co-moving wavenumber vector. Eventually when a(t) oo, all effective DM 
dissipates into effective cosmological constant Aq for any initial conditions of Xi. 

In summary F{K,,\) gravity gives particularly simple equations. The meaning of these 
equations has been explored in part in the VA model of Zhao (2007). The model has the 
effect of a non-uniform Dark Energy fluid, which mimicks galactic Dark Matter, can seed 
cosmic perturbations, and mimics the effect of a coupled Dark Matter and Dark Energy 
fluid. 



4.2. Homogenous and isotropic universe 

As a second check, we consider the general case of F' ^ est, but in the simple case of 
the expanding uniform universe. The only non-zero component of the Einstein tensor are 

• Goo = 3i?^ 

• — Gyy — Gzz — ~(o^ + 2aa) 



We have therefore the 00 modified Einstein equation 



3(1- aF') H'^ + ^M^F = 8nGp (79) 



and the modified pressure equation 



l-)^aF'\- + OiF' E --^^^-kGP (80) 

These results are identical to that of Zlosnik ct al. (2007). This means simply that the C4 
term does not contribute to the expansion except for providing a zero point of pressure. 



4.3. A Possible origin of Cosmological Constant 



To see the C4 term can contribute as cosmological constant, let us consider Hubble 
expansion in the simple case where we set a = 0. For such models /C = 0. The equations for 
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expansion become very simple 

3H^ + ^— = SttGp (81) 
-H'-2^-^ ^ 8nGP, (82) 

Qj Zi 

so the equation of state of the vector field is 

w — —\ for the Hubble expansion at all redshift (83) 

Following Zhao (2007) we set the zero-point F{}Csoiar) = in solar-system like strong 
gravity regime where fCsoiar ~ lO-*^^, since the gravity near the earth's orbit is about lO^M, 
where M ~ 10~^°m/sec^, so 



2Xd)C. (84) 



K. 



olar 



Taylor expand in the limit of weak gravity /C ~ 0, wc have ^^^^ ~ — Aq + "^A/C ~ — Aq -|- 
-^di(j)di(j)^ which has no first order term, but can have a zero point constant Aq, given by 

Ao = -Wcl^ / 2\dK, ~ In K,,oiar, (85) 

*^ ^solar 

where for reasons evidently later we take /i — 1 — X — ^ p^f^- As we will see Aq plays the 
role of the cosmological constant. Interestingly Aq ~ M^lnlO^^ ~ 36M^ ~ Hq, which is 
the observed amplitude of the cosmological constant. The logarithm factor explains why the 
observed A is significantly greater than M^. 



4.4. Static limit 



As another application, we apply our equations to the regime of quasi-static galaxies. 
We set the background expansion factor a — 1. In the static hmit, the spatial terms of the 
vector appear only at second order in all the equations. T^^ has no cross-terms (up to linear 
order), so we find ijj — —(f) from Gij — equation. And the only non-zero component of the 
Einstein tensor is 

-'d^<j> (86) 



Gq — Goo — 2a 



For the vector field we have 



rpO 
^0 



rpx 



a 



2\di. 



rpx 

^X 1 



-FM'^ 
2 



-Ao, 



(87) 
(88) 
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where the pressure term is generally much smaller than Tq ~ SnGp, so the equation of 
state of the vector field is 

w ~ 0, in static galaxies where \k\'^(j) » Aq. (89) 

Prom the Einstein 00th equation, and neglect the pressure term, we find the modified 
Poisson equation 

a-^di {2fxdi(f)) = SttGp, // = 1 - A (90) 
We hence recover eq. (9) of Zlosnik et al. (2007), except that we do not require C4 = 0. 

The above equation resembles the MOND Poisson equation in static limit. However, 
MOND also requires for present-day galaxy — > ^ when y = <^ 1 and /x — > 1 

when y 3> 1, where we identify M with the MOND critical acceleration oq, i.e, y/M = oq ~ 
10~^°m/sec^. With no loss of generality we set C4 = 2, ci = 0. The easiest way to match the 
MOND function with F' together is to require a — 0, y — ]C/2, and 

(91) 

The latter corresponds to the standard /i function of classical MOND, which fits rotation 
curves of hundred nearby spiral galaxies extremely well. 



5. Possible co-variant dependance of the MONDian behavior on redshift, 

environment and history 

As a final application, we note that it is possible to deviate from MOND when wc 
consider galaxy models with a = Ci+ 3c2 + C3 7^ in a non-static universe. As before we set 
H = \]ol+2 • However, the kinetic scalar K, is up to 2nd order 



Hizf „ 1 
'0 



K ~ lOOa-^ + 2y, y= ^ (9,09,0) , (92) 



where /Co = SaH^/M^ ~ lOOa for M ~ Hq/G, Hence we find /x = \l .^t+S^nTi^^iH-M 



l+\y+bQaH{z)^/Hl\ 

depend on redshift. 

Finally, coming back to F{Q) models, the free function now depends on Q = CqM~'^R-\- 
/C, which depends the Ricci scalar, which is crudely speaking the density of the system. Por 
galaxies in an expanding universe, 

Q^cqM i? + /C ~ (6co + 3q;)— + Qo, Qo = , (93) 
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Setting a — —2cq we can also opt out the term or the redshift-dependancy and make 
Q^Qq. E.g., if the MOND function = l-^^g = V^^' MONDian behavior will 
depend on density. The zero gravity 9^0 = region has ~ Q = Qo = n^f^ ~ ~ 
lOOco^p <S 1, where 5p is the over-density over the cosmic mean, and we assume Cq <S 1. So 
the dark matter effect could be bigger in a fluffy galaxy cluster than in a dense galaxy 
in these models. In the solar system Q is big due to high density and strong gravity, hence 
/i = 1, we recover GR like behavior. The F{Q) models contain also a correction to the 
Einstein equation due to a source proportional to —CoWa/B ~ —CoF'Ra/3 ~ 0, where the free 
function F' ~ — 1. This correction can be neglected in the case cq <^ 1 as in most F{R) 
gravity models. 

Coming back our general Lagrangian L{lp, A) with a dynamical (p freedom if = di = 
ds = and d2 = 1. The term A"'A''Va^Vhip creats a quintessence like source term in 
cosmology but does not contribute to static galaxies. However, in time-dependant systems, 
this couphng of A" and cp means that the MOND = in these models has not reached its 



steady state prediction e.g., = y Instead it must be solved from its own equation of 
motion in an unrelaxed system under merging. 

In short, the co- variant version offers new possibilities of tailoring the MOND behavior as 
function of enviornment and redshift and history. These possibilities of covariant dependance 
of the MOND //-function are generally welcome, since some of the MOND's worst outhers 
are with gravitationally lensed galaxy clusters under merging at modest redshift, e.g., the 
Bullet Clusters eA z — 0.3; clusters have generally lower density than spiral galaxies, where 
the empirical formula of MOND applies well. In this sense, the empirical MOND formula is 
not a universal rule, and there are a range of possible fundamental rules giving the effects of 
Dark Matter and Dark Energy. 



6. Conclusion 

We have outlined a framework for studying the dark matter and dark energy effects 
of a vector field. We have isolated a few simple cases where the perturbation equations 
for structure formation are the simplest. Our equations reduce to the non-linear Hubble 
equation and the nonlinear Poisson equations in the literature; Our simplest model with 
Ci 7^ = Ci = C2 = C3 = Cq = di = d2 = d?, = d4 gives particularly simple Einstein's 
Equations. Including other coefficients lead a a range of new behaviors in structure formation. 
We itemize our main results as following. 

• The rotation curves of most spirial galaxies can be explained if we adopt the MOND 
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dielectric parameter //(/C) = 1 — £4_ci^/ _ ^ j^;^^ where K, ~ 2y, where ^ is the 
gravity measured in units of the acceleration scale M. 

• The metric-tracking vector field is described by a four vector Aa = (—1 — 0, Ai, A2, A3). 
It acts as a dark fluid of certain four- velocity. This fluid is able to store up perturbations 
in vacuum in a cold dark matter fashion without being dissipated by photons, hence 
giving the seed for formation of baryonic structures after the epoch of last scattering 
(Dodelson & Ligouri 2006, Zhao 2007). 

• This dark fluid has a non-constant equation of state parameter w. In for the pure C4/C 
case the fluid behaves as aw — —1 cosmological constant Aq in Hubble expansion, and 
w — dark matter in static galaxies. 

• The small amplitude of vacuum pressure Aq ~ Hq is explained by the vector field's 
pressure in galaxies, if the zero point of the pressure is set at the solar system. Here 
Aq is the maximum pressure difference between very strong and very weak gravity. 

• There are co-variant models F{R + }C) models with a = C1 + 3C2 + C3 7^ 0, and/or Cq 7^ 
which allows the MOND dielectric function to depend on redshift, density, hence it is 
no longer a universal rule. 

Our perturbation equations can be fairly straightwardly generalised by superimposing 
two F{Qi) and F{Q2) terms together. E.g., in the VA model (Zhao 2007), one replaces 
F F(c4/C4) -|- F2{J) where F2 ex J oc /C2 in the matter-dominated regime. This J-term 
has effects orthogonal to that of the /C4 term. It can mimic the effects of dark matter in the 
Hubble equation, but does not contribute to galaxy rotation curves. 

The generality of the equations presented here gives the opportunity of exploring var- 
ious realistic cases. With these it is in principle possible to numerically simulate structure 
formation and Cosmic Microwave Background to falsify these F(R + /C) class of models in 
the style of Skordis et al. (2006), and Li et al. (2007). 

This work is part of AH's master thesis project in ENS Paris, done in collaboration 
with HSZ at University of St Andrews. HSZ acknowledges partial support of PPARC Ad- 
vanced Fellowship and National Natural Science Foundation of China (NSFC under grant 
No. 10428308). 
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A. Linear Models with F{IC) = IC 



As a first apphcation of our results, we generahze the hnear model of Lim (2005) to 
include a C4/C4 term. 



We let F(/C) = JC, hence F' = 1. We have 



C4 



+ .. 



The vector equation gives the Lagrange multiplier 



L*{1 + 20) = 3 (ci + C2 + C3) {H' + 2HiP) + 3c2 — + H(t>- - ^dfct) (A2) 



C3 



+2(ci + C2 + c,)^diA, - "-^^dt (a,A) 



The 00 component of the stress-energy tensor 
T — 



„^ _ + 3a//^ + QaHdti^ - i^M^/C + 2^49,^ + ^^^9*(9,A), (A3) 



The spatial diagonal term 



lAn a) 4- -( 



T^^ = -a^^^ -dt{a^a) + Ja^(l + 2ij)M^)C + aa^ -6i/^ + - d^^p (A4) 



6 6a 



The Ox component of the stress tensor 
To, = (C4 - ci) {m<l> + Hdi<f>) + ^^^dt (ai,.) + adt (H) A, + ^^^d, {d,A - 9,>1JA5) 



The spatial off-diagonal terms 



T..^-'-^{H + dt)d^,A,), 



13 



(A6) 



This preprint was prepared with the A AS IM^iX macros v5.2. 
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where the bracket means symmetric permutation of i and j. 

It is reassuring that the above equations agree with those of Lim (2006) if we set C4 = 0. 
This confirms our results up to the hnear order in the case that F' — 1. 



